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We extend the theory of Euler integration from the class of con- 
structive functions to that of "tame" : -valued functions (definable 
with respect to an o-minimal structure). The corresponding integral 
operator has some unusual defects (it is not a linear operator); how- 
ever, it has a compelling Morse-theoretic interpretation. In addition, 
we show that it is an appropriate setting in which to do numerical 
analysis of Euler integrals, with applications to incomplete and un- 
certain data in sensor networks. 

definable functions | o-minimal structures | sensor networks | euler character- 
istic 

Integration with respect to Euler characteristic is a homo- 
morphism J x -dx ■ CF(X) — > Z from the ring of con- 
structive functions CF(X) ("tame" integer- valued functions 
on a topological space A) to the integers Z. It is a topologi- 
cal integration theory which uses as a measure the venerable 
Euler characteristic \- Euler characteristic, suitably defined, 
satisfies the fundamental property of a measure: 

X (AuB) = x(A) + x(B)-x(AnB), [1] 

for A and B "tame" subsets of X . We extend the theory to 
R-valued integrands and demonstrate its utility in managing 
incomplete data in, e.g., sensor networks. 



into simplices. The Mayer- Vietoris principle asserts that x is 
a measure (or 'valuation') on O, as expressed in Eqn. pp. 

The Euler integral is the pushforward of the trivial 
map X i > {pt} to f x d X : CF(X) -» CF({pt}) ^ Z satisfy- 
ing J x 1a dx = x(A) for 1a the characteristic function over a 
definable set A. From the definitions and a telescoping sum 
one easily obtains: 

hd%= J2 s x{h = s} = J2x{h> s}- X {h<s}. 

- C s = -co s = 

[2] 

Because the Euler integral is a pushforward, any definable 
map F : X -> Y induces F, : CF{X) -> CF(Y) satisfying 
J x hdx = j Y F t hdx- Explicitly, 

F*h(y) = / hd X , [3] 
Jf-Hv) 

as a manifestation of the Fubini Theorem. 

The Euler integral has been found to be an elegant and 
useful tool for explaining properties of algebraic curves [4] and 
stratified Morse theory [22, 5 , for reconstructing objects in 
integral geometry [21j . for target counting in sensor networks 
[TJ, and as an intuitive basis for the more general theory of 
motivic integration |8] [7] . 



Constructible integrands 

Because the Euler characteristic is only finitely additive, one 
must continually invoke the word "tame" to ensure that x 
is well-defined. One means by which to do so it via an O- 
MINIMAL structure [23], a sequence O = (C„) of Boolean 
algebras of subsets of R n satisfying a small list of axioms: 
closure under products, closure under projections, and finite 
decompositions in 0\. Elements of O are called definable 
sets and these are "tame" for purposes of integration theory. 
Examples of o-minimal structures include (1) piecewise-linear 
setsQ (2) semi-algebraic sets; and (3) globally subanalytic 
sets. 

Definable functions between spaces are those whose graphs 
are in O. For X and Y definable spaces, let Def (X, Y) de- 
note the class of compactly supported definable functions 
h : X — * Y, and fix as a convention Def(A) = Def(A, R). 
Let CF(X) = Def (X, Z) C Def (X, R) denote the ring of CON- 
STRUCTIBLE functions: compactly supported Z- valued func- 
tions all of whose level sets are definable. Note that in gen- 
eral, definable functions (even definable 'homeomorphisms') 
are not necessarily continuous. 

We briefly recall the theory of Euler integration, estab- 
lished as an integration theory in the constructible setting 
in [151 1201 1211 124] and anticipated by a combinatorial ver- 
sion in [3] 1111 1121 119] . Fix an o-minimal structure O on 
a space X. The geometric Euler characteristic is the func- 
tion x '■ C ~ * Z which takes a definable set A 6 O to 
X(A) = dim iff M (A;R), wne re H? M is the Borel- 

Moore homology (equivalently, singular compactly supported 
cohomology) of A. This also has a combinatorial definition: 
if A is definably homeomorphic to a finite disjoint union of 



(open) simplices \J then x{A) = £,(-1) 



Algebraic 



Real-valued integrands 

We extend the definition of Euler integration to R-valued in- 
tegrands in Def (A) via step-function approximations. 

A Riemann-sum definition. 

Definition 1. Given h £ Def (A), define: 



L 



hYd X \ 



lim — / \nh\dx- 



1 



[4] 
[51 



h \dx\ = lim - / \nh]dx- 

We establish that these limits exist and are well-defined, 
though not equal. 

Lemma 1. Given an affine function h 6 Def(a) on an open 
k-simplex a , 



h [d X \ = (-l) fc inf h; / h \dx] = (~l) k sup ft. [6] 



Proof: For h affine on a, x{[ n ^i > s } ~ (~l) fe f° r all 
s < n infer h, and otherwise. One computes 

If 1 °° 

lim — / \nh\dx = lim — > x{L n ^J > s } — ( — 1) inf ft. 

J " s = 

The analogous computation holds with xlT 72 ' 1 ! > s } = ( — l) fc 
for all s < nsup CT h, and otherwise. ■ 

This integration theory is robust to changes in coordi- 
nates. 



topology asserts that x 1S independent of the decomposition 



1 Some authors require an o-minimal structure to contain algebraic curves, eliminating this particular 
example. 
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Lemma 2. Integration on Def (A) with respect to \_dx\ an d 
\dx] is invariant under the right action of definable bijections 
ofX. 

Proof: This is true for Euler integration on CF(X); thus, it 
holds for J x [nh\ d\ and J x \nh~\ d\. ■ 
Lemma 3. The limits in Definition^ are well-defined. 
Proof: The tri ANGULATION theorem for Def (A) [23] states 
that to any h G Def(X), there is a definable triangulation (a 
definable bijection to a disjoint union of open affine simplices 
in some Euclidean space) on which h is affine on each open 
simplex. The result now follows from Lemmas Q] and [5] B 

Integrals with respect to [dx\ and \dx~\ are related to to- 
tal variation (in the case of compactly supported continuous 
functions). 

Corollary 1. If M is a 1- dimensional manifold and h G Def(M) 

is continuous, then 

[ h [d x \ =- f h \d x ] = itotvar(ft). [7] 

Jm Jm z 

Proof: Apply Lemma [T] to an affine triangulation of h which 
triangulates M with the maxima {p{\ and minima {qj} as 0- 
simplices and the intervals between them as f-simplices. To 
each minimum qj is associated two open f-simplicies, since M 
is a 1-manifold. Thus: 

/ h [d x \ = J2 + " 2 = ^totvar(Ti). 

Jm j j l 

This equals — J M h \d\~\ via an analogous computation. B 

This result generalizes greatly via Morse theory: see 
Corollary [5] One notes that [d\\ and \d\~\ give integrals 
which are conjugate in the following sense. 
Lemma 4. 

J h\d X ] = - J -h[d X l- [8] 

Proof: Apply Lemma [T] to an affine triangulation of h, and 
note that sup CT h = — inf CT — h. B 

The temptation to cancel the negatives must be resisted: 
see Lemma [S] below. 

Computation. Definition [T] has a Riemann-sum flavor which 
extends to certain computational formulae. The following is 
a definable analogue of Eqn. [2] . 
Proposition 2. For h G Def(A), 



/ h[d X ] = / X {h>s}- X {h<-s}ds [9] 

r r oc 

/ h\dx\ = / X {h> s}- X {h<s}ds. [10] 

JX J s=0 



Proof: For h > affine on an open fc-simplex a, 

/ h ld x \ = (-i) k Mh= / x (<r n {h > s})ds, 

and for h < 0, the equation holds with — X (o~ H {h < — s}). 
The result for J \dx\ follows from Lemma |4] B 

ft is not true that J x h [dx\ = J °° sx{h = s}ds: the 
proper Lebesgue generalization of Eqn. [2] is the following: 
Proposition 3. For h G Def (A), 



/ h [dx\ = lim - [ s X {s <h < s + e}ds [11] 
Jx <^o+ 6 Jm 

[ h\d X ] = lim - / sx{s < h < s + e}ds. [12] 
Jx e J R 



Proof: For h affine on an open fc-simplex a, and < e suffi- 
ciently small, J R s x{s < h < s + e} ds = e ( — f ) fc (— | + inf a h) 
and f R sx{s < h < s + e} ds = e (-f) fc (— f + sup CT h). B 



Morse theory. One important indication that the definition of 
J [dx\ is correct for our purposes is the natural relation to 
Morse theory: the integrals with respect to [rfxj an d \dx~\ are 
Morse index weighted sums of critical values of the integrand. 
This is a localization result, reducing from an integral over all 
of X to an integral over an often discrete set of critical points. 

Recall that a C 2 function h : M —* R on a smooth mani- 
fold M is Morse if all critical points of h are nondegenerate, in 
the sense of having a nondegenerate Hessian matrix of second 
partial derivatives. Denote by C(h) the set of critical points of 
h. For each p G C(h), the MORSE index of p, fx(p), is defined 
as the number of negative eigenvalues of the Hessian at p, or, 
equivalently, the dimension of the unstable manifold W u (p) 
of the vector field — V/i at p. 

Stratified Morse theory !f0 is a powerful generalization 
to triangulable spaces, including definable sets with respect 
to an o- minimal structure [5] [22] . We may interpret [dxj and 
\dx\ in terms of the weighted stratified Morse index of the 
graph of the integrand. 

Definition 2. For X C R" definable and h G Def (A), define the 
co-index ofh, T*h to be the stratified Morse index of the graph 
ofh, C A xR, with respect to the projection n : A xR — > R: 



(Th)(x)= lim x\BAx)r\{h<h{x) + e}), [13] 



where B e (x) is the closed ball of radius e about x G A. The 
index T t is the stratified Morse index with respect to height 
function —iv : i.e., X#h = h) or 



(l,h)(x)= lim x (B e (x) n {h > h(x) - e'}) . [14] 

e'<§Ce->0+ V / 

Note that 1,,T : Def (A) -> CF(X), and the restriction 
of these operators to CF(X) is the identity (every point of a 
constructible function is a critical point). The two types of 
integration on Def (A) correspond to the Morse indices of the 
graph with respect to the two orientations of the graph axis 
— the projections it and — n. 

Theorem 4. For any continuous h G Def(X), 

f h[d x \= [_hThd x ; / h\dx\= [_hl»hd x . 
Jx Jx Jx Jx 

[15] 

Proof: On an open fc-simplex a C A C R n in an affine tri- 
angulation of h, the co-index T*h equals (— l) dlm (°") times the 
characteristic function of the closed face of a determined by 
inf CT h. Since h is continuous, J_hT*hdx ~ {— f) dnn ' <T ^ inf CT h. 
Lemma [1] and additivity complete the proof; the analogous 
proof holds for X* and \dx~\ ■ B 

Corollary 5. If h is a Morse function on a closed n-manifold 
M , then: 
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Proof: For p £ C(h) a nondegenerate critical point on an n- 
manifold, Th(p) = and l,h(p) = ■ 

From this, one sees clearly that the relationship between 
[dxj an d \d X ~\ is regulated by Poincare duality. For exam- 
ple, on continuous definable integrands over an n-dimensional 
manifold M, 

[ h [dxl =(-!)" / h[d X \. [18] 

JM JM 

The generalization from continuous to general definable 
integrands is simple, but requires weighting T*h by h directly. 
To compute f x h \d\\ , one integrates the weighted co-index 

lim h(x + e')x ( BAx) D {h < h(x) + e'\) [19] 

e'<CE-»0+ V / 

with respect to d\. 

Corollary [5] can also be proved directly using classical 
handle-addition techniques or in terms of the Morse complex, 
using the fact that the restriction of the integrand to each 
unstable manifold of each critical point is unimodal with a 
unique maximum at the critical point. It is also possible to 
express the stratified Morse index — and thus the integral 
here considered — in terms of integration against a charac- 
teristic cycle, cf. [101 153] . 

One final means of illustrating Corollary [5] is to use a de- 
formation argument. Let h be smooth on A and <j> t be the flow 
of —Vh. Then the integral is invariant under the action of <f>t 
on h; yet the limiting function hoo = rimt_>oc> ho(j> t is constant 
on stable manifolds of — X7h with values equal to the critical 
values of h. We have not shown that the limiting function is 
constructible (this depends on the existence of definable in- 
variant manifolds — we are unaware of relevant results in the 
literature) and thus do not rely on this method for proof but 
rather illumination. 

The integral operator 

We consider properties of the integral operator(s) on Def (A). 

Linearity. One is tempted to apply all the standard construc- 
tions of sheaf theory (as in [201 EH) to J x : Def(A) -» R. 
However, our formulation of the integral on Def(X) has a 
glaring disadvantage. 

Lemma 5. J x : Def (A) — > R (via \_dx\ or \d X ~\) * s n °t a ho- 

momorphism for dim A > 0. 

Proof: 

1=/ l[d X \^ [ x[dx\+[ (l-x)[d X \=l + l = 2. 

J [0,1] J [0,1] J [0,1] 

■ 

This loss of functoriality can be seen as due to the fact that 
[f + g\ agrees with [/J + \_g\ only up to a set of Lebesgue 
measure zero, not X" meas ure zero. The nonlinear nature of 
the integral is also clear from Eqn. [15] , as Morse data is 
non-additive. 

The Fubini Theorem. In one sense, the change of variables 
formula trivializes (Lemma [2j . The more general change of 
variables formula encapsulated in the Fubini theorem does 
not, however, hold for non-constructible integrands. 
Corollary 6. The Fubini theorem fails on Def (X) in general. 

Proof: Let F : X = Y\\Y ~ * Y be the projection map 
with fibers {p} U{p}. Any h 6 Def (X) is expressible as / TJ g 
for /, g £ Def(y). The Fubini theorem applied to F is equiv- 
alent to the statement 

/ / +/ 9= f h= f F*h = [ f + g 
Jy Jy Jx Jy Jy 



(where the integration is with respect to [dxj or \dx~\ as de- 
sired). Lemma [5] completes the proof. ■ 

Fubini holds when the map respects fibers. 

Theorem 7. For h e Def(X), let F : X -> Y be defin- 
able and h-preserving (h is constant on fibers of F). Then 
J Y F, h [d X \ =J x h [d X \ , and J Y F, h \d X ] = f x h \d X ] . 

Proof: An application of the o-minimal Hardt theorem [23] 
implies that Y has a partition into definable sets Y a such that 
F~ 1 (Y a ) is definably homeomorphic to U a x Y a for U a defin- 
able, and that F : U x Y a — ► Y a acts via projection. Since h 
is constant on fibers of F, one computes 

/ F*h[d X \= f h X (U a )ldxl= [ Kd X \. 

JY a JY a JU a xY a 

The same holds for f \dx] ■ B 

Corollary 8. For h G Def(X), J x h = J R h*h. In other words, 

f h[d X \ = / sx{h = s}[d X \, [20] 
Jx Jm 

and likewise for \dx \ ■ 



Continuity. Though the integral operator is not linear on 
Def(X), it does retain some nice properties. All properties 
below stated for J [dx\ hold for J \dx~\ via duality. 

Lemma 6. The integral J [dxj : Def(X) — » R is positively ho- 
mogeneous. 

Proof: For / £ Def(X) and A G R + , the change of variables 
variables s As in Eqn. [9j gives J Xf [dx\ = A J f [dx\ ■ B 

Integration is not continuous on Def(X) with respect to 
the C° topology. An arbitrarily large change in J ^[^xj ma y 
be effected by small changes to h on a (large) finite point set. 
In some situations the "complexity" of the definable functions 
can be controlled in a way sufficient to ensure continuity. 

One example arises in the semialgebraic category. Fix a 
(finite) semialgebraic stratification S of a compact definable 
X, and consider definable semialgebraic functions algebraic 
with respect to this stratification (that is such that the re- 
striction of the function to any stratum S £ S is a polynomial 
Ps). The resulting linear space (filtered by the subspaces 
of polynomials of bounded degree) can be equipped with the 
structure of a Banach space, by completing the family of semi- 
norms ||P||s,k = maxsgs ||Ps||c™ , where n = dimX. Then 
fx ' L^xJ becomes a continuous (non- linear) functional on this 
Banach space. The proof results, essentially, from the Bezout 
theorem (mimicking Thom-Milnor theory): the total number 
of critical points graph of a polynomial of degree Dona fixed 
semi-algebraic set is bounded by 0(D n ). The generalization 
to increasing (refined) stratifications is straightforward. 

Integration itself defines a natural topology for Def (X) on 
which integration is continuous. Define the L 1 e-neighborhood 
of h G Def (A) as the intersection of the C° e-neighborhood 
(definable functions with e-close graphs) with those functions 
g £ Def (A) satisfying I J" / — g [dx\ \ < e. This provides a 
basis for an L 1 topology on Def (A). As a consequence of 
Lemma [J] the definition is independent of the use of [dxj or 
\dx\- 

The interested reader may speculate on other function 
space topologies on Def (A). 
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Duality and links. There is an integral transform on CF{X) 
that is the analogue of Poincare-Verdier duality [22] ■ It ex- 
tends seamlessly to integrals on Def(X) by means of the fol- 
lowing definition. 

Definition 3. The duality OPERATOR is the integral transform 
V : CF(X) -» CF(X) given by 

Vh(x) = lim / hl Bt{x) d X , [21] 
e^0+ J x 

where B ( is an open metric ball of radius t. 

We extend the definition to V : Def(X) -> Def (X) by in- 
tegrating with respect to \d\\ or \d\\, interchangeably, via: 
Lemma 7. T>h is well- defined on Def(X) and independent of 
whether the integration in i21]) is with respect to [d\\ or \d% \ ■ 
Proof: The limit is always well-defined thanks to the Conic 
Theorem in o- minimal geometry [23]. To show that it is inde- 
pendent of the upper- or lower-semicontinuous approximation, 
take e > sufficiently small. Note that by triangulation, h 
can be assumed to be piecewise-affine on open simplices. Pick 
a point x in the support of h and let {a} be the set of open 
simplices whose closures contain x. Then for each i, the limit 
hi(x) := \im y -> x h(y) for y G ct; exists. Then, applying Eqn. 
[21 j . one computes 

Vh{x) =Y J {-l) dim ° i hi{x), [22] 

i 

independent of the measure [dxj or T^xl • ' 

For a continuous definable function h on a manifold M, 
T>h = (— 1) m h, as one can verify by combining Eqns. [5] 
and [21] . This is commensurate with the result of Schapira 
[20] that V is an involution on CF(X). 
Theorem 9. Duality is involutive on Def(X): T> o T>h = h. 
Proof: Given h, fix a triangulation on which h is affine on 
open simplices. From Eqn. [22] , we see that the dual of h at 
x is completely determined by the trivialization of h at x. Let 
L x h be the constructible function on B t {x) which takes on the 
value hi(x) on strata cr; C\B € (x). (Though this is not necessar- 
ily an integer-valued function, its range is discrete and there- 
fore it is constructible.) As L x h is close to h in B t (x) (this 
follows from the continuity of h on each of the strata), T>h is 
close to DL x h in B e (x): indeed, the total Betti number of in- 
tersections of strata with any ball B e (y) is bounded, and Euler 
integral of a function small in absolute value is small as well. 
Hence the definable function T> 2 h is close to the constructible 
function V x h with e small. As T> 2 L x h(x) — L x h(x) = h(x), 
the result follows. ■ 

One can define related integral transforms. For example, 
the LINK of h G CF{X) is defined as 

Ah(x) = lim / hl dBtix) d X . [23] 

6^0+ J X 

The link of a continuous function on an n-manifold M is mul- 
tiplication by 1 + ( — 1)™, as a simple computation shows. In 
general, A = Id — T>, where Id is the identity operator. 

Integral transforms 

Integration with respect to Euler characteristic over CF(X) 
has a well-defined and well-studied class of integral trans- 
forms, expressed beautifully in Schapira's work on inversion 
formulae for the generalized Radon transform in d\ [2T]. In- 
tegral transforms with respect to \dx\ an d \dx~] ar e similarly 
appealing, with applications to signal processing as a pri- 
mary motivation. Examples of interesting definable kernels 



for integral transforms over Euclidean R" include \\x — y\\, 
(x,y), and g(x — y) for some g. These evoke Bessel (Han- 
kel) transforms, Fourier transforms, and convolution with g 
respectively. The choice between [dxj an d \dx~] makes a dif- 
ference, of course, but can be amalgamated. Example: for 
fixed kernel K, one can consider the mixed integral transform 
h i — ► J x hK [dx\ — J x hK \dx~] ■ In the case of K(x, f ) = (x, £), 
this transform takes 1a for A compact and convex to the 
'width' of A projected to the £-axis. 

Convolution. On a vector space V (or Lie group, more gener- 
ally), a convolution operator with respect to Euler character- 
istic is straightforward. Given f,g G CF(V), one defines 

(/*<?)(*)= / f(t)g(x-t)d x . [24] 
Jv 

Convolution behaves as expected in C'F(V). By re- 
versing the order of integration, one has immediately that 
J v f * g dx — J v fdxJ v gdx- There is a close relationship 
between convolution and the Minkowski sum, as observed in, 
e.g., [TT]: for A and B convex and closed 1a * 1b = 1a+b, 
cf. [241 120) . Convolution is a commutative, associative opera- 
tor providing CF(V) with the structure of an (interesting [J|) 
algebra. 

Convolution is well-defined on Def (V) by integrating with 
respect to [o^xj or T^xl • However, the product formula for 
f f * g fails in general, since one relies on the Fubini theorem 
to prove it in CF(V). 

Linearity. The nonlinearily of the integration operator pre- 
vents most straightforward applications of inversion formulae 
a la Schapira. Fix a kernel K 6 Def(A x Y) and consider 
the integral transform Tk '■ Def(X) — * Def(Y) of the form 
(Ti<;h)(y) = J x h(x)K(x,y)[dx\(x). In general, this operator 
is non-linear, via Lemma [5] However, some vestige of (posi- 
tive) linearity survives within CF. 

Lemma 8. The integral transform Tk is positive-linear over 
CF+{X) = Def(A,N). 

Proof: Any h G CF + (X) is of the form h = J2 k a ^u k for 
a k G N and U k G Def(A). For h = 1 A , T K h = J A Kld X \- Ad- 
ditivity of the integral in [dx\ (via Eqn. [T]) combined with 
Lemma [6] completes the proof. ■ 

This implies in particular that when one convolves a func- 
tion h G CF + (R n ) with a smoothing kernel (e.g., a Gaussian) 
as a means of filtering noise or taking an average of neighbor- 
ing data points, that convolution may be analyzed one step 
at a time (decomposing h). 

Integral transforms are not linear over all of CF(X), since 
f —h[dxi 7^ — J h[dx\- However, integral transforms which 
combine [dx\ and \dx~\ compensate for this behavior. Define 
the measure [dx] to be the average of [rfxj an d T^xl : 

Jx h[dx] = \(l x h[dX] + Jx HldX] ) ' [25] 
Theorem 10. Any integral transform of the form 

(T K h)(y)= f h(x)K(x,y)[d X ](x) [26] 
Jx 

is a linear operator CF(X) — > Def(V). 

Proof: From Lemma [8] T is positive-linear over CF + (X). 
Full linearity follows from the observation that J x —h[dx] = 
— f x h[dx], which follows from Lemma [1] by triangulating h. 
■ 

As a simple example, consider the transform with kernel 
K{ x >£) = i x iO- The transform of 1a with respect to [dx] 
for A compact and convex equals a 'centroid' of A along the 
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£-axis: the average of the maximal and minimal values of £ 
on OA. Note how the dependence on critical values of the 
integrand on dA reflects the Morse-theoretic interpretation of 
the integral in this case. 

Integration with respect to [dx] seems suitable only for 
integral transforms over CF. On a continuous integrand, the 
integral with respect to [dx] either returns zero (cf. the in- 
tegral of Rota Q3|]) or else the integral with respect to L^xJi 
depending on the parity of the dimX, via Eqn. [18] . 




Fig. 1. Shown is a sampling of an upper semicontinuous constructible integrand 
h : — > N, sampled over a random network with holes. Values of the integrand 
are color-coded by height, white being zero. One third of the sample points sample h 
with an error of ±1, uniformly distributed. The Euler integral of this sampling with 
respect to d\ is 64, a poor approximation to the true J h d\ = 9. 




Fig. 2. The above sampled integrand, averaged over nearest neighbors. The in- 
tegral of this IR-valued integrand with respect to [^xj ' s approximately 9.52: a 
reasonable approximation to the true J h d\ = 9. 



Applications of definable Euler integration 

The Euler calculus on is quite useful; the extension to Def 
deepens this utility and opens new potential applications, of 
which we highlight a few. 

Sensor networks. An application of Euler integration over 
CF(X) to sensor networks problems was initiated in [T|. Con- 
sider a space X whose points represent target-counting sensors 
that scan a workspace W . Target detection is encoded in a 
SENSING relation S C W x X where (w,x) € S iff a target 
at w is detected by a sensor at x. Assume that sensors count 
the number of sensed targets, but do not locate or identify 
the targets. The sensor network therefore induces a TARGET 
counting function h : X — > N of the form h = Y2 a lt/„j 
where U a is the target support — the set of sensors which 
detect target a. Euler integration allows for simple enumera- 
tion: 

Theorem 11. (Q]) Assume h G CF(X) and x{U a ) = N ^ for 
all a. Then the number of targets in W is precisely i J x hd\- 

Since the target count is presented as an integral, it is pos- 
sible to accurately estimate the answer when the integrand h 
is known not on all of X (a continuum of sensors being an 
idealization) but rather on a sufficiently dense grid of sample 
points (physical sensors in a network). 

The R-valued theory aids in establishing expected values 
of target counts in the presence of confidence measures on 
sensor readings. Let M = {xi} denote a discrete set of sensor 
nodes in R n , and assume each sensor returns a target count 
h(xi) £ N and a fluctuation measure c(xi) € [0, 1] obtained, 
say, by stability of the reading over a time average. View h 
as a sampling over Af of the true target count / = ~}2 a ljy Q . 
Assume that nodes with fluctuation reading have perfect 
information (h = / at xt) and that c correlates with error 
|/ — h\. Assume that sensor nodes Af are the vertex set of a 
triangulation T. 

The integral of an extension of / over a triangulation gives 
a terrible approximation to J hd\'- an error of ±1 on if nodes 
can cause a change in the integral of order K. More specif- 
ically, if h = f + e, where e : J\f —> { — 1,0,1} is an error 
function that is nonzero on a sparse subset M' C M, then, for 
certain infelicitous choices of Af', \fh — ff\ = |A/"'|. 

A R-valued relaxation can mitigate errors by using fluctu- 
ation c as a weight. Let N(i) be a collection of neighboring 
nodes to Xi, where neighborhood can be defined via distance 
(if available) or edge-distance (in an ad hoc network or trian- 
gulation). Define h to be the result of averaging the value at 
Xi 6 TV over N(i), with c as a weight. Specifically, 

~ E wgJ y(i) c(y)h(y) 

2^y£N{i) ( '\y) 

This nearest-neighbor convolution damps out local variations. 
The resulting integral with respect to [dx\ will tend to mit- 
igate localized errors, thanks to the Morse-theoretic formula: 
see Figs. 1-2. More numerical investigation is warranted. 

Such averaging leads to non integer-valued integrands. By 
using integration with respect to [dx\ or \dx] for upper/lower 
semi-continuous integrands associated to an averaged signal h, 
one obtains an expected value of J hdx- This can be particu- 
larly illuminating when a network has incomplete information, 
e.g., a hole. Holes in a network can be modeled by setting the 
confidence measure c to zero and averaging. 

Statistics and mode counting. The previous applications lend 
themselves to more general statistical ends. Consider a 
smooth distribution / : 1" — > [0, oo) of compact support 
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and bounded variation. The problem of mode-counting — of 
decomposing / into a convex combination of unimodal sum- 
mands — is of interest to statisticians, even in the univariate 
case [9] [2]. Gaussian summands are commonly used, though 
not exclusively [131 114] . The techniques of this paper are rel- 
evant (see [18] for Morse structures associated to mixtures of 
multivariate normal distributions). 

One way to interpret mode-counting is as a topological de- 
convolution problem. Assume that / is of the form f — u * h 
for h £ CF + (R n ) and u a unimodal distribution supported 
on a neighborhood of € R n , where unimodal means that all 
upper excursion sets {u > c} are contractible for c > 0. One 
assumes that h is of the form h = X^iLi ^-u t where the supports 
Ui are sufficiently simple so that 11*1^ has no self-interference 
(the support of u * strongly deformation retracts to Ui). 

The mode-counting problem is equivalent to computing 
J hdx given / and some information about u, say, its height 
maxw = J u[dx]- If the convolution formula for CF held over 
Def, then the number of modes would be J /L^xJ divided by 
Umax- However, the non-linearity of the integral with respect 
to Yd\\ precludes this solution. Indeed, this nonlinearity mir- 
rors the interaction of unimodal summands in a distribution. 
Just as two modes can interfere, creating a local maximum 
when an increasing and a decreasing portion of the modes 
are summed, the Euler integral over Def loses linearity when 
increasing and decreasing integrands are combined. 

The development of good algorithms for integral trans- 
forms over \d\\ or even [dx] will be useful not only for mode- 
counting, but also for explicit mode decomposition, since the 
extraction of the modes themselves from the unknown factor 
h involves deconvolving / = u * h. 

Numerical integration. Though integration with respect to 
Euler characteristic has a lengthy history, there appears to 
be no treatment of numerical integration, even in the sim- 
pler setting of CF + (M rl ). The central problem (in the con- 
structive and definable categories) is how to estimate f x h 
given the values of ft on a discrete subset of X. As in the case 
of numerical integration for Riemann integrals, one typically 
assumes something about the features of h and/or the density 
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and extent of the sampling set. In [I], the present authors 
give a formula for estimating J hdx given a discrete sampling 
of h £ CF(R 2 ) which correctly samples connectivity data of 
excursion sets. This formula generalizes to the definable cat- 
egory: 

Proposition 12. For h £ Def (R 2 ) continuous, J h [dxi = 
/ Po{h >s} + (3 {h > -s} - f3 {h < s} - (3 {h < -s} ds, 

J 3 = 

[28] 

where /3o(-) = dim Ho(-; R), the zeroth Betti number. 
Proof: Apply the homological definition of x to Eqn. [9]; then, 
use Alexander duality in the plane to reduce all terms to /3o 
quantities. ■ 

The value of Proposition [12] is that it allows for computa- 
tion based on /3o quantities. Such connectivity data are easily 
obtained from a discrete sampling via clustering. We have 
implemented this formula in software (see Fig. ??). However, 
for more general integration domains than R 2 , duality formu- 
lae are less helpful. One general result on refinement follows 
from continuity of the integral operator. 

Theorem 13. For h £ Def(A) continuous, let hpL be the 

piecewise-linear function obtained from sampling h on the ver- 
tex set of a triangulation T of X . As the sampling and trian- 
gulation are refined, 

lim / h PL [d X i = / h[d X ]. [29] 

This result relies crucially on continuity and does not ap- 
ply to CF(X). A more desirable result would be a measure 
of how far a given sampling is from the true integral. This 
seems challenging. We note that the Morse-theoretic formu- 
lae [16]- [17] allow one to reduce the domain of an integral to 
a (typically finite) set of critical points. This 'focusing' prop- 
erty of integration over Def should be a starting point for good 
numerical algorithms, especially for integral transforms. 
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